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The CDF collaboration has set lower limits on the masses of the Z' bosons occurring in a 
range of E§ GUT based models. We revisit their analysis and extend it to certain other 
E§ scenarios as well as to some general classes of models satisfying the anomaly cancellation 
conditions, which are not included in the CDF analysis. We also suggest a Bayesian statistical 
method for finding exclusion limits on the Z' mass, which allows one to explore a wide range 
of the U(l)' gauge coupling parameter. This method also takes into account the effects of 
interference between the Z' and the SM gauge bosons. 



1 Introduction 

A neutral Z 1 gauge boson appears in numerous models containing the SM gauge symmetry group 
along with an additional U(l) symmetry (for a review, see^). Gran d Unified Theories (GUTs) 
larger than the original SU(5) model, such as 50(10) ^ or £? 6 1214| ; b rea k down to the SM as: 
E 6 ->• 50(10) x 17(1)^, with in turn 50(10) -> 5*7(5) x U(l) x ^ SM x U(l)fi The Z' 
thus surviving at the electroweak (EW) scale can be written as the linear combination, 

Z' = cos a cos (3Z X + sin a cos (3Zy + sin f5Z$. (1) 

If kinetic mixing^ *^^ with the hypercharge group U(1)y is neglected by setting a = in 
eq. ([I]), one obtains some well-known Z' bosons by adjusting (3. Th ese include Z x {j3 = 0°), Z$ 
{p = 90°), Z v rj -52.2°)EI, Zi (p ps 37.8° P, Z s (fi ~ 23.3° jMUi an d z N (/3 ps 75.5°)H21. The 
inclusion of kinetic mixing results in certain other phenomenologically interesting cases, such as 
Z dp h (a, P ~ -78.5°, 37.8°) which does not couple to the d-type quarks, Z R (a, P « 50.8°,0°)E1 
which couples to the right-handed fermions only and Zb-l (a, P ~ —39.2°, 0°)^^, where B is the 
baryon number and L the lepton number of an ordinary fermion. The Zlr which exists in models 
with left-right symmetry^^ is equivalent to the linear combination: \/3/5 (a Zr — Zb-l/2&), 



where a = J g\l g\ cot 2 6w — 1, with 6\y being the weak mixing angle and gL,R being the 
SU(2)L t R coupling strengths, respectively. The Zlr studied here corresponds to the specific 
case of = gn- In addition to these Eq based models, a Z s t r ing from a specific superstring 
modelQl 

and a sequential Zs m are also included in this analysis. 

A number of other classes of one or more-parameter models have been discussed EUSED. 
Without the assumption of unification at the GUT scale, but assuming nullification of anomalies 
using three families of exotics (which is not the case in some super symmetric models), there are 



Table 1: fermion charges, e/, in terms of Zf, their values in the corresponding ^-parameter models. Also given 

are the angles a and /3 yielding these models. 
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four 'one-parameter' models^', with fermion charges B — xL, d — xu, q + xu and 10 + x5, with 
x arbitrary. The last three of these correspond to the generalized Eq charges in eq. ([1]). (q + xu 
corresponds to arbitrary superpositions of Y and B — L, while 10 + x5 corresponds to the Eq 
models without kinetic mixing.) We have, therefore, normalized the fermions charges, Zf, in the 
x-parameter models to their Eq values, ej. These charges are given in Table [TJ 

2 Z' production at CDF and limits on its mass 

The total cross-section for the Drell-Yan (DY) process at a hadron collider, with a neutral 
gauge boson B as the mediator and fi + fi~ as the outgoing particles, is given as^ 

a = - [ * MdMa diff , (2) 
s Jo 

where M is the invariant mass of the muon pair, y/s is the center-of-mass (CM) energy, and 
/ dx 1 -Y / HM 2 )^\ff(x 1 ,M 2 )f?(x 2 ,M 2 ) + f((x 1 ,M 2 )f^(x 2 ,M 2 )}, (3) 

JM 2 /s x l „ iVc 

where A^ c = 3 is the quark color factor. x\ and x 2 (= above are the momentum fractions 

of the ingoing partons having parton distribution functions (PDFs) j^!^ for hadrons A and B. 
a s is the strong coupling constant and K = Kq + Ke, with Kq being the QCD iC-factor and 
Ke being the multiplicative factor due to QED corrections. a{M 2 ) is equal to 
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(\A LL \ 2 + \A RR \ 2 ) (1 + cos^ f + {\A LR \ 2 + \A RL \ 2 ) (1 - cos 0*) 2 j , (4) 



with 9* being the polar angle defined in the CM frame and the individual amplitudes given as 

A " ~ "* 2 + Mi-^liM z r z C >^<» + M>-MlliM z ,T z , C ?^ C ?^ < 5 » 

where i,j are run over L, R. Q and e are the electric charges of the contributing quark and the 

z z 1 z z 1 

muon, respectively. C L ' R (/) = g\^j^ R (/), with g\ = ej sin 8 W cos 6 W being the gauge coupling 

strength of the Z boson, g2 the Z' coupling and e^' R the EW and U(l)' charges of the fermion 
/• Tz z' are t ne total decay widths of the Z and Z' bosons having masses Mz z'- 

We first follow the CDF analysis!^ an d use th e lo expression given in eq. ([2]) to calculate 
the cross-section due to the Z' alone, neglecting the 7 and Z contributions to the amplitudes in 
eq. ©. We employ LO CTEQ6L PDFs^and mass-dependent NNLO i^c^and NLO K E ^ 
values, and assume that the Z' decays into SM fermions only, which are taken to be massless. 
We achieve up to 99.8% agreement on the 95% confidence level (C.L.) Mz> lower limits with 
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Figure 1: Contours in Mz> limits mentioned on/beside them, x shows the location of a particular Eg model 
within a contour. The dotted, dashed and dot-dashed lines correspond to the three ^-parameter models. Z v * is 
a leptophobic boson which will not be observable in the di-lepton channels. 



Table 2: Numerical limits in GeV on the mass of the Z' boson in various models. 
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Table 3: Limits on M z > in the Z x model from the LHC. 
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the CDF analysis, for the E§ models included therein, and obtain new limits for the rest of the 
models. The numerical values of the limits are given in Tableland the a, (3 parametrization of 
the models with contours in Mz> is plotted in Fig. [TJ In Table [3] we give limits obtained for the 
LHC with a similar approach for some expected integrated luminosity and CM energy values. 

3 Bayesian statistical method 

The CDF analysis uses signal templates generated with a fixed resonance pole width, T = 
2.8% x Mz>- However, there is no fundamental reason to only look for such a narrow Z' . A 
wide Z' resonance, implying a strongly coupling boson, could well be scattered over a few bins 
and no significant enhancement above the background will be visible. Besides, the effects of 
interference between the various bosonic contributions to the propagator, i.e., between 7, Z and 
Z' (see eq. ©), are lost in their approach. These effects could in principle cause a considerable 
enhancement or dip in the number of events in several accompanying low M _1 bins, e.g., in the 
case of a strongly interacting Z' boson with mass just beyond the kinematic reach of the CDF. 
Finally, the CDF limits assume a fixed GUT-based 52 and it is not straightforward to extend 
the limits to other g<i values, particularly in the strong coupling regime (see, however^J) . 

Therefore, we propose a Bayesian statistical method which allows one to vary 52 in order to 



obtain the corresponding limit on Mz>- It is based on the likelihood function L, written as 



Hm = nf P(ni\m), (6) 

where P is the Poisson probability of finding nj events given fa expected events in the ith bin 
with B total bins. L in eq. ([6]) is then evaluated using two hypotheses: the null hypothesis, 
L(/i b |n), assumes that the DY process occurs only via 7 and Z, and the signal hypothesis, 
L(/Lt*|n), with /i* = fj, b + [/,', wherein the Z' boson also contributes to the cross-section along with 
the SM gauge bosons. The SM events, expected in an invariant mass bin are calculated as 

2EC f 

fi b = Cctsm = / AdMp 1 / MdM^{M- l \M- x )(T d if f + n BG , (7) 

s Jbin JO 

where C = 2.3 fb _1 is the integrated luminosity at the CDF, E = 0.982 is the detector efficiency 
and A is the CDF acceptance, which is a mass-dependent multiplicative factor. M p in the above 
equation is the muon-pair mass measured by the detector, u-qg refers to the non-DY events and 
the probability density is given as 

p(M- 1 |M- 1 ) = M p b a e- b /T(a), (8) 

with a = (M^/A) 2 , b = JW~ 1 M~ 1 /A 2 , where A = 0.17 TeV" 1 is the variance. The purpose of 
the above probability function is to smear over the DY background before distributing it into 
bins, hence accounting for the mis-identification of an event in a bin where it does not actually 
belong. For fjf, ax 2 function is constructed as 

4, = -2LLR = -2 ln(M^l ) = _ 2 ^ _ ^ + ni ^(M), ( 9) 

and is minimized to obtain the best-fit values in the given range of g% and Mz' , which correspond 
to a Z' boson with cross-section // best favored by the data. The contours in #2 and M^, 
for a certain C.L. value specified by the allowed number of standard deviations, A^ 2 , from 
the minimum, can then also be drawn, giving the exclusion limits on these parameters. Our 
preliminary contours are given in Fig. [2]for Z x as a representative model. The numerical value of 
the 95% C.L. limit is 913 GeV for g 2 = 0.461, which is about 21 GeV higher than the CDF value. 
We next plan to undertake a global analysis including constraints from electroweak precision 
dataEl Eventually, a similar statistical analysis of the LHC data, as soon as it is released, will 
also be performed. 
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Figure 2: Exclusion contours for the Z x boson. The dotted, dot-dashed and double-dot-dashed curves correspond 
to the C.L. values given and the solid curve represents the CDF limits generalized to other values. 
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